CCTP-2011-24 



F-theorem, duality and SUSY breaking in 
one-adjoint Chern-Simons-Matter theories 



Takeshi Morita and Vasilis Niarchos 

Crete Centre for Theoretical Physics, 
Department of Physics, University of Crete, 71003, Greece 

takeshi@physics.uoc.gr, niarchos@physics.uoc.gr 
Abstract 

We extend previous work on M = 2 Chern-Simons theories coupled to a single adjoint chi- 
ral superfield using localization techniques and the F-maximization principle. We provide 
tests of a series of proposed 3D Seiberg dualities and a new class of tests of the conjec- 
tured F-theorem. In addition, a proposal is made for a modification of the F-maximization 
principle that takes into account the effects of decoupling fields. Finally, we formulate and 
provide evidence for a new general non-perturbative constraint on spontaneous supersym- 
metry breaking in three dimensions based on Q-deformed S 3 partition functions. An explicit 
illustration based on the known analytic solution of the Chern-Simons matrix model is pre- 
sented. 
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1. Introduction 

Important progress has been recently achieved in understanding some of the key prop- 
erties of three-dimensional supersymmetric conformal field theories (SCFTs). The main 
technical tool has been localization which (after proper regularization) allows to express the 
free energy of the Euclidean CFT on the three-sphere 

F = -\og\Z s ,\ (1.1) 

in terms of a compact matrix integral expression [TH3]. 

The free energy F provides a variety of useful information. For example, three-dimensional 
SCFTs with M = 2 supersymmetry have a conserved U(l) R-symmetry that sits in the same 
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supermultiplet as the stress-energy tensor. In a general interacting SCFT this symmetry re- 
ceives quantum corrections and the quantum numbers associated with it become non-trivial 
functions of the parameters that define the theory. The proper coupling of the matter fields 
to the curvature of S 3 makes F a function of the R-charges. It has been proposed that the 
exact U(l) R-symmetry of the theory is that R-symmetry which locally maximizes F [2J. 
Non-trivial tests of this proposal have appeared in [H-TTO]. 

Secondly, there are by now many known examples of three-dimensional SCFTs that ex- 
hibit dualities analogous to the S-duality and Seiberg duality in four dimensions. Matching 
the free energy F of dual theories provides non-trivial checks of these dualities [TTHTo] (analo- 
gous checks of 3D dualities based on the computation of superconformal indices, see e.g. [T5] . 
have been discussed in pTflTg]). 

Moreover, it has been noted that there is a close relation between superconformal indices 
in four dimensions [TUll2"0] and sphere partition functions in three dimensions [2"TH2"3"] . In 
this context, the mathematical identities that relate the superconformal indices of dual 
theories in four dimensions are closely related to the mathematical identities that relate 
the sphere partition functions of dual theories in three dimensions. This observation raises 
the possibility of a promising 3D-4D connection. 

Finally, it has been proposed [HJI21] that F is a quantity that plays in three dimensions 
the same role that the c-function plays in two dimensions ~\ In analogy to Zamolodchikov's c- 



theorem in two dimensions [27] (or the g-theorem for one-dimensional boundary RG flows 
|2"5]). an F-theorem has been conjectured in three dimensions, which states that F decreases 
along RG flows. Successful checks of this proposal have appeared in [TUlEniEI] • 

In this paper we will discuss all three of the above aspects in the context of a special 
class of M = 2 Chern-Simons-Matter (CSM) theories defined as Af = 2 Chern-Simons 
theory at level k coupled to a single chiral superfield in the adjoint representation of the 
gauge group, which we take to be U(N). We consider this theory in the presence or absence 
of superpotential interactions. This is an interesting theory for the following reasons: 

(1) In the absence of superpotential deformations the theory is believed to be exactly 
superconformal for all N, k [32J and exhibits an R-symmetry with large quantum cor- 
rections. At strong coupling, e.g. large 't Hooft coupling A = N/k 3> 1 in the large- N 
limit, the R-charge of the adjoint superfield tends to zero [S3]. As a result, with in- 
creasing coupling more and more fields hit the unitarity bound and decouple from the 
rest of the theory as free fields, thus truncating the interacting part of the chiral ring 



lr The F- theorem was first argued as a general aj-theorem in the context of holography in 
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from below. 

(2) In the presence of superpotential deformations and sufficiently strong coupling the 
theory exhibits spontaneous supersymmetry breaking [34]. In the range of parameters 
where a supersymmetric vacuum exists the superpotential deformation gives rise in the 
infrared (IR) to a new class of interacting fixed points. These theories are believed to 
exhibit a Seiberg-like duality that acts on the theory in a self-dual manner [31]. Having 
a minimal matter content these examples provide some of the simplest illustrations 
of Seiberg duality in three dimensions and thus a useful playground for attempts to 
understand the general underpinnings of such dualities in field theory. 

Our main goal will be to analyze the implications of these features on the properties of 
the free energy F. A preliminary analysis of F-maximization in these theories in the large- 
N 't Hooft limit was presented in [EJIH]- Evaluating F in the saddle point approximation 
it was confirmed numerically that F-maximization provides at A <C N (but not necessarily 
A <C 1) results consistent with a set of non-perturbative inequalities derived in [S3]- This is a 
non-trivial test of the validity of the F-maximization principle. In section [2] we supplement 
this result with a more thorough examination of the solutions of the saddle point equations. 
Besides the one-cut solutions of Refs. [ElE] we demonstrate the existence of a large class 
of multi-cut solutions and find numerical evidence that the one-cut solutions of Refs. [5],[H] 
are in fact the dominant ones (namely, they have the smallest free energy in the large-iV 
limit). The multi-cut solutions are expected to play an interesting role in the analysis of the 
non-perturbative (in 1/N) instanton contributions to the matrix model that expresses the 
free energy after localization. 

At A ~ N Ref. [S] observed numerically that one of the non-perturbative bounds of [33] 
gets violated by a naive application of the F-maximization principle. A potential source for 
this violation is the sizable number of free decoupled fields in this regime. In section 0] we 
test this expectation by suitably modifying the F-maximization principle to account for the 
decoupling fields. We find that the modification leads to an R-symmetry consistent with 
the non-perturbative bounds of [33] • We propose that the implemented modification is the 
general way to deal with the accidental symmetries of decoupling fields in F-maximization. 

Another missing link in the F-maximization prescription is the following. Ref. [2] argues 
for F-extremization, but the more precise statement about maximization (namely the sign 
of the second derivative of F with respect to the trial R-charges) is currently mostly an 
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empirical observation rather than a generally argued facto In the one-adjoint theories that 
we consider we always observe numerically a single maximum that reproduces correctly in 
appropriate regimes independent information about the exact R-symmetry. 

In addition, in section [2] we observe regions in (A, i?)-space (A being the 't Hooft coupling 
and R the trial R-charge of the adjoint chiral superfield), where the one-cut solutions cease 
to exist. This effect occurs precisely when the supersymmetric vacuum of the physical 
(superpotential-deformed) theory disappears. In the special case of the Chern-Simons matrix 
model, that expresses the free energy of the topological M = 2 Chern-Simons theory, we 
demonstrate that this effect is directly related to the vanishing of the full non-perturbative 
partition function. 

This observation motivates the possibility of a general relation between zeros of the Q- 
deformed sphere partition function, which leads to localization, and spontaneous supersym- 
metry breaking. We formulate and further motivate this potential relation in very general 
terms in section [3] and appendix [Bl As a full-fledged theorem this relation would provide a new 
non-perturbative constraint on spontaneous supersymmetry breaking in three-dimensional 
QFTs that supplements the known constraints arising from the Witten index. 

In section we consider the second theme outlined in the beginning, namely Seiberg 
duality. In a simple non-topological case, that involves a U(l) and a U(2) CSM theory, we 
are able to compute the free energies analytically and demonstrate that they match exactly 
in the dual theories. In more general cases, we find convincing numerical evidence for this 
matching in the large-iV limit using our saddle point solution^ In subsection 15.31 further 
aspects of the duality of interest are discussed in some length. 

Finally, the one-adjoint CSM theories possess a large class of RG flows induced by su- 
perpotential interactions. In section we outline this web of RG flows and the general 
predictions implied by the conjectured F-theorem. In subsection 16.31 we find numerical ev- 
idence for the validity of these predictions in a set of examples with increasing complexity, 
thus providing new tests of the conjectured F-theorem. 

In two extra appendices (app. [X]and[C]) we summarize for the benefit of the reader a set 
of useful results that are alluded to, but not explicitly analyzed, in the main text. 

2 Besides the explicit verification of maximization in specific examples, a general relation between F- 

maximization and volume-minimization has been proposed for large- N CSM gauge theories with candidate 

Sasaki-Einstein gravity duals [1HB] . See also the more recent work [33] . 

Preliminary signs of potential disagreement with Seiberg duality reported in [5] are not verified. 
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3 


4 


5 


6 


7 


8 


9 


10 


A n+1 





0.35 


0.65 


0.9 


1.15 


1.4 


1.7 
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Table 1: The numerically determined values of \* n+l using F -maximization [B] forn = 3, 4, . . . , 10. 
Summary of notation 

The one-adjoint CSM theories that we consider are characterized by two integers, N the 
rank of the U(N) gauge group, and k the level of the CS interaction^ Following [33] we 
frequently denote the theories without any superpotential interactions as A. New interacting 
fixed points arise by deforming the action with the superpotential interactions W n+ i = 
TrX n+1 , n = 1,2, . . .. We will variably denote the resulting theories either as A n+ i or in 
more detail as U(N)^ +1 \ 

A large- iV 't Hooft-like limit is possible and defined as 



N, k ->■ oo , A 



N 
~k 



fixed . 



1.2) 



In Refs. j8J[33] it was argued that as we increase the coupling A in the undeformed theory A 
the R-charge function R(X) decreases monotonically towards zero attaining the values 

2 



R(X* n 



n+U 



n 



1 



n 



1,2, 



;i.3) 



at the critical couplings A* +1 . At each of these points the operator TrX n+1 is, by definition, 
marginal. Using F-maximization, [BlEj, one finds for the first few values of n the values of 
A* +1 listed in Table [Q above. 

It has been conjectured that the A n+1 theories exhibit a Seiberg-like duality [34]. In 
the RG flows that give rise to the A n+ i theories the deforming operator Tr X n+l is relevant 
in both the 'electric' and 'magnetic' versions of the theory when A G [A* +1 ,n — A* 



n+U- 



Following the nomenclature of analogous situations in four- dimensional gauge theories, we 
will sometimes call this interval the 'conformal window 7 of the A n+1 theory. A review of the 
precise details of Seiberg duality in the A n+ i theories appears in section 



2. S 3 partition function in the large-A limit 

The main object of interest in this paper is the S 3 partition function Z s z. Using local- 
ization techniques [IH3] one can re-express Z53 as a matrix integral. In our case, this matrix 



We assume k > 0. The case of negative k can be obtained by a simple parity transformation. 
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integral reads 

i . / N \ N N 

Z s3 (N,k,R) = ± [lle^Utj n(2sinh(7r^)) 2 n e^' R+ ^ = e -WM> (2.1) 
where ij are matrix eigenvalues running over the real line and tij = U — tj. £(z) is the function 




The free energy that we maximize with respect to R is 

F=\{F + F). (2.3) 

In the A n+1 theories the R-charge is fixed by the superpotential interaction W n+ i = 
Tr X n+1 , so in order to compute the free energy F in these theories one simply sets R = ^-j- 
in the expression (12. ip . 



2.1. Saddle point equations and their solutions revisited 

In the large- N limit the main contribution to Z#3 comes from saddle point configurations 
that obey the system of algebraic equations 



J^l , (1 - R) smh(27r^) + Uj sin(27nR) 

CO til yTrt j 



Xi = -u-\ — y 

1 A 1 N^. 



cosh(27rtjj) — cos(27ri?) 



0, z = l,2,...,iV. 

(2.4) 



At a saddle point configuration 



N - at N N 



F(X, N) = - log TV! + £ ±^t 2 + £ log (4 sinh 2 (7rf;;)) + £ ^(1-R+ iUj) ■ (2.5) 

i=l i<j i,j=l 

In general, the N tj's that solve these equations are complex numbers (here C- valued func- 
tions of the parameters R, A). 

In lack of an efficient analytic method of solving these equations Refs. [HJE] proceeded to 
analyze them numerically. In Ref. [5] solutions were found by introducing a fictitious time 
coordinate r (following J3H]) and considering the dynamical evolution described by the set 
of differential equations 

a d ^ = X l . (2.6) 

GST 

In these equations a = e ine (8 G [0, 2ir)) is a free constant whose choice affects (in some cases) 
the solution that one converges to. Another crucial input that affects the obtained solution 
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Figure 1: Plot (a): A typical one-cut distribution of the eigenvalues ti (i = 1, ...,N) in the 
complex plane. This particular plot has been obtained for N = 100, A = 1 and R = 0.28. Plot (b): 
A typical two-cut solution depicted here for N = 100, A = 1.5 and R = 0.5. The specific example 
crosses the imaginary axis at ±|. 



are the initial conditions at r = 0. Using two independently written codes in Mathematica 
and Fortran we searched for different solutions varying a and the initial conditions. 

Besides the one-cut solutions that were reported in Refs. [HUH] we have been able to detect 
many more multi-cut solutions with branches crossing the imaginary axis through the points 
_l_mi^ j_ (R+m)i ^ ( m = ^ 2, . . .). A typical one-cut solution appears in Fig. [TYa) and a typical 
two-cut solution in Fig. W(p). The multi-cut solutions can also be traced analytically in the 
weak coupling regime (A<1) with a perturbative computation (some of the key details of 
this computation are summarized in appendix [AT) . 

The one-cut solutions were the ones used for F-maximization in Refs. [S1IS], hence it is 
important to know if they are also the dominant ones in the large- N limit (i.e. if they are 
the ones having the minimal free energy for fixed A, R). We confirmed this fact numerically 
in all cases that we checked. Further credence in the dominance of the one-cut solutions 
is provided by the following facts: (i) at weak coupling they are the ones that reproduce 
the perturbative field theory computation of the R-charge [9], (ii) at finite coupling they 
reproduce to a good accuracy the predictions of Seiberg duality and F-theorem (see below). 

The multi-cut solutions are expected to play an interesting role beyond the saddle point 
approximation. In general, multi-cut solutions are closely related to the non-perturbative 
instanton corrections that control the large-order behavior of the 1/N expansion pJTHlU] . 
Since such effects are beyond the immediate scope of the present work in what follows we 
will concentrate mostly on results based on the one-cut solutions. 
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2.2. One-cut solutions and the functional behavior of F(X, R = in the A n+ i theories 

Now we come to a more interesting observation. We focus on the A n+1 theories and 
fix the R-charge to the value dictated by the deforming superpotential, hence R = 
We would like to understand the functional behavior of the localized free energy function 
F(X, ^Ti) as we vary the 't Hooft coupling A. In particular, we would like to understand 
what happens to this function when A > n, i.e. when the supersymmetric vacuum is lifted. 
In that case the Q-deformed S 3 partition function computed via localization is no longer 
equal to the physical partition function of the A n+ i theory, but we can consider it anyway 
(for example, it gives information about the partition function Zg3(X, R) of the undeformed 
A theory before F-maximization at R = ^-j- — see below for more comments related to this 
aspect). 

A particularly instructive example that is amenable to full analytic control is provided 
by the A 2 theory, i.e. the case with n — 1. In that case, the chiral superfield X is massive 
and the theory flows in the IR to the topological N = 2 Chern-Simons theory. The sphere 
partition function (12.11) simplifies considerably and can be computed exactly for any N, k. 
Setting R = 1 in fl2.ll) we obtain (using the identity i(—z) = —£(z), z G C) 

|Z 5 3(A, J R = l)| = ^ y 

(2.7) 

For the second equality we used a result that is explained, for example, in appendix B of 
Ref . pQ . Eq. (12. 7p gives the absolute value of the partition function of a well-studied matrix 
model (the U(N) Chern-Simons matrix model [JT|) 

1 r ( N 2 \ N 
Z ™ = wJ ( IK 2 ^"^ J n ^smh^,)) 2 . (2.8) 

In the conventional case S is a real parameter. In our case 5* = 2-7riA. 

We would like to draw attention to the following property of the expression (12.71) : the 
exact function Z S 3(X,R = 1) vanishes identically in the supersymmetry-breaking regime, 
Accordingly, in this regime the free energy F, derived by localization, diverges. 
Since our information about the general n case is currently limited to the leading-order 
large-iV limit, it is worth elaborating further on how this divergence exhibits itself in the 
saddle point approximation. 

5 This fact was also observed in Ref. [T2] . 




n ^dt 3 



N 



n(2sinh(vrty))' 

i<j 



N-l 



f,N/2 



m=l 



[ 2 sin 



Tim 



N-m 



k J 



S 



One way to derive the free energy F at leading order in 1/N in the 't Hooft limit is by 
directly taking the large- iV limit of the exact expression (12 .7p . By converting a sum to an 
integral we obtain 

f = ~]og\Z sa \ = ^ 1 d//(//-l)log(2sin(7rA^)) . (2.9) 
For A < 1 this integral evaluates to 

1 = ~ 3A + 2A2) " ik? ( Ll3 ( e "* iA ) " c(3) ) ■ (2 - 10) 

Although not immediately obvious, the rhs of this equation is a real number. For A > 1, i.e. 
in the SUSY-breaking regime, part of the integrand in (12. 9|) lies on a logarithmic branch cut 
and the expression becomes ill-defined. 

Notice that / vanishes at the SUSY boundary value A = 1. Hence, the crossover from 
a well-defined regular / at A < 1 to a diverging / at A > 1 happens discontinuously (more 
comments on the nature of this discontinuity below). The vanishing of / at the general SUSY 
boundary values N = k is a property of the exact \Z$3\ for any N, as can be seen from eq. 
(12. 7p . This generic feature in the A n+ i theory is a natural property given the vanishing of 
the rank of the gauge group in the Seiberg dual theory. 

Another derivation of eq. (12.91) can be obtained by solving directly the saddle point 
equations. An analytic solution of these equations is known (321II2]- For real S in (I2.8P the 
solution is a one-cut solution along the real axis with the eigenvalues occupying the interval 

C = [-c,c] , c = ^arccosh (e s/2 ) . (2.11) 

The resolvent of the solution is 



7T 



u(t) = 2 log 
and the eigenvalue density 



e ae 71 " 4 (cosh(7rt) - 1\J e s - cosh 2 (7rt) 



(2.12) 



1 1 cosh(7rt) — i\/e s — cosh (irt) 
p{t) = --lmu(t) = -— log V 1 I . (2.13) 

71 711 cosh(7rt) + iye 5 - cosh 2 (7rf) 

By analytically continuing this result to S — > 27riA we obtain the one-cut solution that we 
are interested in. From the endpoint of the cut at 

c = -arccosh (> iA ) (2.14) 

it is apparent that as we increase A from to 1, the cut rotates away from the real axis 
towards the imaginary axis. At A = 1 the cut lies fully along the imaginary axis. Including 
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the complex conjugate configurations that correspond to A G [—1, 0] one covers the full range 
of solutions provided by the analytically continued version of (I2.12p . In other words, there 
are no one-cut solutions that can be identified in this way for |A| > 1, which corresponds to 
the non-supersymmetric regime of the A 2 theory. 

This picture is fully reproduced by the numerical results of the previous subsection. We 
reproduce the specifics of the eigenvalue distribution, e.g. the A-dependence of the one-cut 
endpoint in ( I2.14p . and the functional dependence of the free energy dictated by the expres- 
sion f)2.9p for A G [0, 1] (see Fig. [2](a))j^ When we try to extend the numerical computation 
above A = 1 we fail to find any regular one-cut solutions and are instead able to obtain only 
multi-cut solutions. In view of the exact result (12. 7p . it is natural to associate the breakdown 
of the one-cut solutions at A > 1 to the non-perturbative divergence of F. 

From a large-iV matrix model point of view one anticipates the following picture. In 
different regions of moduli space the matrix model is dominated by different saddle points. 
In cases where a one-cut phase dominates the corrections coming from multi-cut phases are 
exponentially small. In other cases, however, different multi-cut phases are of the same 
order, they sum in a coherent way and the partition function can have zeroes. This type 
of asymptotics and phase transitions were studied in Ref. [H|. A more thorough discussion 
on the vanishing locus of the partition function can be found in the unpublished lecture 
notes 05] E 

Having a detailed understanding of the n = 1 case, we can now proceed to ask what 
happens at generic n where our knowledge is unfortunately restricted only to a numerical 
solution of the large- N saddle point equations. 

Our data indicate that the picture is qualitatively similar at all values of n. Nice one- 
cut solutions exist throughout the whole supersymmetric interval, A G [0,n], but disappear 
disco ntinuously when the supersymmetric vacuum of the A n+ i theory is lifted. Extending the 
lessons of the n = 1 case on the basis of the general picture outlined above, we propose that 
this breakdown is always due to the non-perturbative divergence of F in the supersymmetry- 
breaking region. 



6 One can view the good comparison of these data with the analytically expected result (|2.9j) as an 

additional test of the accuracy of our numerical results at N = 100. 

7 In [45j it is pointed out that zeroes of the partition function appear in many statistical systems as 

generalization of Lee- Yang zeroes. By the general theory of asymptotic expansions, they condense along 

anti-Stokes lines and in matrix models they are related to the zeroes of the ^-function (see, for example, the 

footnote in page 39 of Ref. [IS]). We are grateful to Marcos Marino for explaining these aspects to us and 

for drawing our attention to the Refs. 
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Figure 2: Plot (a) depicts the numerically determined renormalized free energy /(A) for the 
topological A 2 theory (blue data points) in the supersymmetric interval [0, 1]. The blue curve depicts 
the analytic result (|2,9|) . Plot (b) provides the analogous numerical result for the A4 theory (blue 
data points) for A G [0, 3]. The dashed red curve is a fit provided by the expression (I2.15|) , 

A specific example of the numerically determined A-dependence of the free energy / in 
the supersymmetric interval is depicted in Fig. W(b) for the A 4 theory There is no known 
exact expression for the localized free energy in this non-topological case. Interestingly, we 
observe that the function 

f f u,n=3 = 0.8 £dn{n - 1) log ^2 sin (^j^j (2.15) 

provides a rather good fit to the numerical data. Qualitatively similar results can be obtained 
also for other values of n with a fit of the form 

ffit,n = c n dfi(fi - 1) log ^2 sin f-^-j 1 , (2.16) 

for an appropriately chosen n-dependent number c n . An appealing feature of this functional 
form is that it obeys the correct Seiberg duality relation (15.131) that will be discussed later in 
section Unfortunately, despite this property, (12.161) cannot be the exact expression of the 
matrix model free energy at leading order in 1/N for generic n. At small A it can be checked 
that /(A, R) has the expansion detailed in eq. (12.171) below, which would imply according to 
( I2.16P c n = 1 for all n. This would be inconsistent with the presented numerical results. 

It remains an interesting open problem to determine analytically the precise expression 
of / at any n and also to verify non-perturbatively in the matrix model that Z s z = when 
N > nk. As a related point, in Ref. [12] it was shown analytically in a U(N C ) Chern-Simons- 
Matter theory with Nf pairs of (anti)-fundamental chiral multiplets that the matrix model 
Z S 3 is always vanishing in the supersymmetry-breaking regime N c > \k\ + Nf. 
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Figure 3: A plot of the free energy of a free field f(R) = -£(l - R). 
We conclude this section with a couple of brief comments. 

The first comment has to do with the observed discontinuity of the free energy F across 
the SUSY-breaking boundary Assuming F (derived by localization) has to diverge in the 
SUSY-breaking region, the discontinuous jump from a finite value (in our set of examples 
F — 0) to +oo is consistent with the intuition that the sphere free energy is a natural measure 
of the number of degrees of freedom. With increasing coupling gauge theories typically 
exhibit a reduction of the number of degrees of freedom, which is indeed in qualitative 
agreement with the observation that F is a monotonically decreasing function of A. In 
that sense, it would have been puzzling to observe that the free energy increases in the 
supersymmetric regime near A = n to meet the diverging behavior at A > n continuously^] 

The second comment is related to the general behavior of the function F(X,R). In 
the context of F-extremization, it is interesting to know whether F has a single or multiple 
maxima as a function of R at fixed A. Numerically we always find a single smooth maximum. 
More specifically, at small fixed A < 0.2 we observe that F(R) exhibits a behavior similar to 
that of the free energy of a free field. Within perturbation theory one finds 

/(A, R) = ±F(\, R) = - l - log(47r 2 A) - 1(1 - R) + 0(X 2 ) . (2.17) 

The contribution £(1 — R) is the free energy of a free field depicted in Fig. [3j As we increase 
A the behavior changes to one where the one-cut solutions cease to exist above a critical 
value of *g That criticaJ va!ue appears to be set by ^«(A) = This is constat with 
the previous discussion about the (non) existence of one-cut solutions in the A n+ i theories. 

8 We thank Andrei Parnachev for a question that prompted us to think about this aspect. 

9 We have not been able to determine numerically the precise manner in which this change takes place as 

we vary A. We consistently observed a critical value of R for all the values of A > 0.5 that we checked. 
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3. Spontaneous SUSY breaking and zeros of the localized S 3 partition function 



The observations of the previous section suggest an intriguing relation between the Q- 
deformed S 3 partition function that leads to localization and spontaneous supersymmetry 
breaking, similar in some respects to the one provided by the Witten index [16] . Here we will 
attempt to formulate and motivate such a relation in a broad descriptive manner postponing 
a more rigorous treatment to future work. 

Consider a general three-dimensional classically superconformal quantum field theory]^ 
Place the Euclidean version of this theory on a three-sphere S 3 and compute the partition 
function 

Z s3 =Je- s - t f^ (3.1) 

where S is the action of the theory, Q is one of the supercharges and V some fermionic 
functional. We assume S is a properly defined action that leads to a correct path integral 
treatment of the theory even at strong coupling, t is the coupling of the Q-exact deformation. 
V should be chosen in such a way that the deformation {Q, V} is positive definite and the 
standard localization technique can proceed without obstructions [17]. For example, in the 
case of general Yang-Mills or Chern-Simons gauge theories with arbitrary matter content 
one can use the prescription of [IH5]. 

In a theory with a supersymmetric vacuum a standard argument shows that the t- 
derivative of Zg3(t) is identically zero. Hence, the physical S 3 partition function defined 
at t = and its infinitely Q-deformed version at t = oo are identical 

Z s3 (0) = hm Z s z (t) = Zf a c) . (3.2) 

In cases where the standard localization techniques [IrE] (based on a UV Lagrangian for- 
mulation) can proceed unobstructed and capture correctly the quantum IR physics, the rhs 
of this equation, Z^° c \ can be further expressed in terms of a matrix integral formula. For 
general M = 2 supersymmetric theories Zg° ^ is a function not only of the parameters of the 
theory (e.g. the 't Hooft coupling etc.), but also of a set of trial R-charges. F-maximization 
postulates that the physical values of the R-charges can be determined by minimizing \Zg° \. 
So far, this is the familiar story that has led to the recent progress in three-dimensional su- 
persymmetric quantum field theories. 

\Zg° |, computed with this specific procedure, is the quantity one would naturally com- 
pute given a specific theory without any a priori knowledge about the fate of supersymmetry 



10 For reasons that will become obvious in a moment we assume that the theory has at least Af = 2 
supersymmetry. 
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at the quantum level. In a theory with spontaneous supersymmetry breaking the t-derivative 
of Z s z{t) is non-zero in general and the equality (13. 2p is no longer valid (see appendix |B1 for 
more details). For that reason the quantity Zg° c ^ does not compute in this case the phys- 
ical S 3 partition function of the theory Zgs(0). Nevertheless, we would still like to claim 
that Zg° QS) is an interesting quantity to consider even when supersymmetry is spontaneously 
broken. 

Our interest in this quantity stems from the potential validity of the following state- 
ments in a general classically superconformal quantum field theory (with at least M = 2 
supersymmetry) : 

(a) Zg° c ^ = implies that the theory exhibits spontaneous supersymmetry breaking. 

(b) Spontaneous supersymmetry breaking implies Zg° = 0. 

Here we refer to spontaneous supersymmetry breaking on S 3 . We believe this is ultimately 
related to spontaneous supersymmetry breaking on other three-manifolds, e.g. S 2 x IR which 
is directly related to the information carried by the Witten index [15] . We will not, however, 
address this connection in this work. Furthermore, since we refer to spontaneous supersym- 
metry breaking in a finite volume we should mention that this does not necessarily imply 
that supersymmetry is also broken at infinite volume and the infinite volume limit may need 
to be treated with care [4"o] . 

The first statement is natural for the following reasons. If the theory possesses a super- 
symmetric vacuum, Zg° c ^ computes, as argued above, the physical partition function on S 3 . 
Then, assuming the free energy F is a finite quantity for a quantum field theory with a finite 
number of degrees of freedom we conclude that |^° c ^| is a strictly positive quantity. The 
assumption that F is finite is natural and expected to hold if F defines a sensible analog of 
the c-function in three dimensions]^ Therefore, the equation Zg° c ^ = cannot be true in 
a theory with a supersymmetric vacuum and we should infer that, when it vanishes, Z^° c ^ 
does not compute the physical partition function Zs^O). In that sense, the conclusion of 
spontaneous supersymmetry breaking seems to be a natural one. 

A potential pitfall should be noted here. Currently, the most efficient method to compute 
Zg° c ^ proceeds via the localization technique which allows to recast the infinitely Q-deformed 
path integral (13.11) in terms of a matrix integral. As we pointed out above, this method relies 



n Iii particular, the application of the .F-theorem would be subtle in RG flows between theories with infinite 
F. In such cases Fjjv = Fir = oo and the validity of the ^-theorem, which states Fjjv > Fir, would be in 
general unclear. 
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on the assumption that a certain UV Lagrangian description of the theory captures correctly 
the quantum IR physics. Incorrect conclusions may be reached if this extra assumption 
fails. For instance, one could encounter situations where a supersymmetric vacuum exists 
but the matrix integral computation gives Z^° c ^ = 0. An example of this sort has been 
presented recently in independent work by Benini, Closset and Cremonesi [IB]- The theory 
in question is three-dimensional super-QCD with gauge group U (N c ) and Nf = N c — 1 pairs 
of (anti) fundamental multiplets. This theory has a quantum-mechanically deformed moduli 
space of supersymmetric vacua, but the naive computation of Z^" ^ via the standard matrix 
integral gives Z^° c ^ = 0. Given the fact that the IR theory is a free theory of neutral chiral 
multiplets, one can argue that the actual value of Z^"^ is non-zero and equal to the value 
of the undeformed partition function in accordance with our general conjecture. The main 
lesson is that the sphere partition function Z$3 should be computed always, and at any 
deformation parameter t including t = and t = 00, using the proper formulation of the 
path integral. 

The reverse statement (b) could be motivated (or eventually argued) in a manner outlined 
in appendix [HI Under a set of assumptions outlined in appendix [B] we argue that there are no 
states that can contribute non-trivially to the localized path integral Zg° c ^ when the theory 
exhibits spontaneous breaking of supersymmetry. This argument is qualitatively similar to 
the one showing that spontaneous supersymmetry breaking requires that the Witten index 
vanishes. As in the case of the Witten index, here also we would have been unable to reverse 
this argument to conclusively argue that Zg° c ^ = implies spontaneous supersymmetry 
breaking (namely statement (a)). It is interesting that, unlike the Witten index case, for 
Z^t we can find an independent motivation for (a) that goes through the finiteness of F 
as a rough measure of degrees of freedom and its relation to the F-theorem. 

If correct, the combined statements (a) and (b) provide a new powerful non-perturbative 
constraint on spontaneous supersymmetry breaking in three dimensions. Under the assump- 
tion that the standard localization techniques [TH3] proceed unobstructed, they postulate 
that the free energy of a matrix model determines whether a corresponding 3D quantum 
field theory does or does not exhibit spontaneous supersymmetry breaking. The examples 
studied in this paper (and the examples in [12]) provide some evidence in favor of such 
constraints. The study of Z^° c ^ in more examples with known patterns of spontaneous su- 
persymmetry breaking may provide further insights into the fate of these statements. We 
hope to return to this task and a more formal argument in favor of (a) and (b) in a different 
publication. 
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Finally, we may add the following potentially interesting aspect to the above story. In the 
context of statements (a) and (b) it may also be interesting to ask whether Zg° c ^ combined 
with the 3D-4D connection, provides a way to relate supersymmetry breaking patterns in 
three and four dimensions. Similarities in such patterns are observed, for example, between 
the four-dimensional M = 1 SQCD theory and three-dimensional M = 2 CSM theories with 
fundamental and antifundamental matter, the four-dimensional M = 1 one-adjoint SQCD 
theories and N = 2 CSM theories of the type considered in this paper with additional matter 
in the (anti)fundamental representations etc. [53] . 

4. Decoupling fields and modifications of the F- maximization principle 

It can be argued [33] that the exact U(l) R-charge R of the adjoint chiral superfield X in 
the A theory tends to zero at strong coupling. For example, in the large 't Hooft coupling 
limit limA->oo F(A) = 0. In [ElE] it was indeed shown with F-maximization that R is a 
monotonically decreasing function of A that seems to asymptote to zero. 

As the R-charge decreases towards zero the scaling dimensions A n+1 = (n + 1)R of the 
chiral ring operators Tr X n+l also decrease and, accordingly, with increasing coupling more 
and more chiral ring operators hit the unitarity bound and decouple as free fields. It is easy 
to show [33] that the general operator Tr X n+l hits the unitarity bound for the first time 
when the operator Tr A" 4(n+1 ) becomes marginal. When a field decouples, a new accidental 
symmetry occurs and the exact U(l) R-symmetry of the theory can mix with it. In general, 
this affects the validity and formulation of the F-maximization principle. 

In the large- iV, finite A <C N computations of Refs. [HUH] these effects were not important, 
because only a small finite number of free fields existed which had a negligible contribution 
to the total free energy of the theory that scales like N 2 . However, as one increases A and 
makes it of the same order as N the increased number of decoupled fields can have a sizable 
contribution to the free energy and the F-maximization principle must either be dropped 
or appropriately modified. The computation of these effects requires a direct analysis of the 
saddle point equations in the 'M-theory' limit — N — > oo, k finite — which is unfortunately 
not a straightforward exercise. In the language of Ref. [B] this difficulty is due to the non- 
cancelation of long-range forces on the eigenvalues. On a more superficial level, Ref. [5] 
observed numerically that the R-charge obtained by naive F-maximization at A ~ N = 100 
violated a non-perturbative bound on the R-charge that follows from Seiberg duality and 
the assumption that R is a monotonically decreasing function of A. A potential source of 
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Figure 4: The numerically determined R-charge curve in the very strong coupling regime for 10 < 
A < 100 and N = 100. The upper blue points represent the result of the standard F -maximization 
recipe. The lower black square points represent the result of the modified F -maximization recipe 
(|4.2|) which is designed to take into account the effects of the decoupling fields. The exact result 
is required by consistency to lie within the bounds of the lower dashed black curve 2 (a+i) an ^ ^ e 
upper dotted red curve 2 \+\ l° wer bound holds for A £ N and the upper bound for A £ N/2) [SJ. 

this discrepancy is the considerable number of decoupling fields. 

The general question of interest here is the following: is there a suitable modification of 
the F- maximization principle that takes properly into account the effects of decoupling fields 
in a theory with large anomalous dimensions? In corresponding situations in four dimensions 
one is instructed to maximize a modified a-function where the 't Hooft anomalies associated 
to the free fields have been subtracted [iniED]. In this section we want to consider and test 
an analogous modification for F-maximization in three dimensions. 

A natural course of action would simply be to subtract the free energy of the decoupling 
fields. If O is a decoupled operator with trial scaling dimension A its contribution to the 
free energy is 

F free = -£(1 - A) . (4.1) 

Therefore, if m chiral operators, Tr X, Tr X 2 , . . . , Tr X m have decoupled in our theory at 
't Hooft coupling A, we propose that one should not maximize the free energy F(X, R) that 
follows from the matrix integral (12.11) . but the modified free energy 

m 

F mod (A, R- m) = F(X, R) + J2 «(1 " mR) . (4.2) 

i=l 

We implemented numerically this modified F-maximization principle to the A theory 
and the result is plotted (for N = 100) in Fig. HI At each A we determined m, the number 
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of free fields, in a self-consistent manner by requiring Tr X m to be the highest in polynomial 
power chiral ring operator that is decoupled after F mo d-maximization is implemented. In 
Fig. H] the blue round dots are the result of the standard .F-maximization procedure that 
exhibits the above-mentioned violation near A ~ 100 (i.e., the curve traced by the blue data 
points crosses the upper-bounding red curve given by ) 0. The black square dots are the 
result of F mo d-maximization. We observe that the modified principle brings down the value 
of R towards the black dashed lower-bound curve. The new R-charge curve lies between the 
upper and lower bound curves and exhibits no signs of bound violationo Moreover, a naive 
fit of the black data points gives at large A the following asymptotic behavior of R(X) 

R(X) -0.46A- - 91 (4.3) 

which agrees with the intuition that R(X) should trace closely (or asymptote to) the lower 
bounding curve given by 2 (i+\) 0- 

This short computation is suggestive of the validity of F mo d-maximization, but unfortu- 
nately does not provide conclusive evidence. The above numerical computation has been 
performed at a finite, but large, value of N where the leading order saddle point approxima- 
tion is corrected by several sources oil/N contributions. It is not completely clear how these 
1/N effects compare with the contributions due to the decoupling fields. For that reason, it 
would be desirable to implement and test our proposal in more examples of CSM theories 
that exhibit decoupling fields, e.g. examples that include extra fields in (anti) fundamental 
representations [83] • In such cases, the contribution of decoupling meson-like operators would 
be of the same order as the leading order free energy. 



5. Duality in A n+ i theories 

In Ref. [33] evidence was provided for a Seiberg-like duality in A n+1 theories] 13 ! At finite 
N, k the duality relates the U (N) A n+1 theory at level k to the U(nk — N) theory at level 
— k; in convenient notation 

f/(iV)i n+1) U(nk-N) { ^ 1] . (5.1) 



12 For a clear explanation of the origin and meaning of the bound curves we refer the reader to [51155], 
13 Rcf. [33] considered generalizations of this duality in theories with one or two adjoint chiral superficlds 
and additional chiral superfields in the (anti)fundamental representations. In Ref. |51j a duality was proposed 
for an SU (2) Af = 2 CS theory at level one coupled to a single chiral superfield in the adjoint and no super- 
potential interactions. A Scibcrg duality in three dimensional Chcrn-Simons-Matter theories was originally 
proposed by Giveon and Kutasov in |52| for a theory with matter in (anti)fundamcntal representations. 
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In this case Seiberg duality acts self-dually and exchanges only the rank of the gauge group 
and the sign of the CS level. In the large- iV 't Hooft limit the duality acts by exchanging 
A -H> n — A (where here we define A = f^r) - I n terms of the number of matter fields involved 
this is one of the simplest known dualities in three dimensions. 

Let us briefly review the arguments in favor of this duality. Originally, this duality was 
motivated in Ref. [M] with a standard D-brane argument based on a type IIB string theory 
setup that involves N D3-branes suspended between n NS5-branes and a single (l,k) 5-brane 
bound state suitably oriented to preserve the right amount of supersymmetry. In order to 
argue for duality one moves the fivebranes through each other and considers the effects of 
this motion on the D3-branes (which is the place where the Chern-Simons-Matter theory of 
interest lives) . To the degree that this brane motion is inconsequential for the IR physics, this 
is an argument in favor of duality in the IR field theory on the D3-branes. A potentially subtle 
point in this procedure has to do with the fact that, because of the particular orientation 
of the branes, the fivebranes have to meet and cross each other in spacetime at a point 
where singular behavior can arise triggering a transition. Such phenomena are known to 
occur, but the presence of the Chern-Simons interaction is believed to alleviate the singular 
behavior and help the theory avoid a transition (see [TJ] for a relevant discussion and explicit 
examples) . 

In our case, additional evidence for the validity of the duality is provided by the following 
facts: 

(1) In the special case where n — 1 the theory flows in the IR to the topological M = 2 
CS theory. Duality in this case reduces to level-rank duality and thus one can check 
explicitly the matching of the S" 3 partition functions [12] . 

(2) In the general n case, we can further deform the superpotential to 

n a 

W = Y - TrX^ 1 "* (5.2) 

and flow in the deep IR to a product of A n . theories (rij < n). Suitably arranging the 
coefficients $ one can flow to a product of A 2 theories 

n 

U(Nt +1) fW)i 2) x • ■ ■ x U(N n )W f £ tf. = JV . (5.3) 

i=l 

(2) 

As already noted, the duality holds in each of the U(Ni) k factors as level-rank duality. 
One might expect that the duality continues to hold as we tune back all the coefficients 
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gi (i — 1, . . . , n) to zero, keeping go non-zero, in order to recover the original U (N)^ +1 ^ 
theory. A more detailed discussion of these RG flows will appear in the next section 
EJ The behavior of the duality under the general superpotential deformation is further 
discussed in appendix O 

In what follows we attempt to find more detailed tests of the duality using the localized 
S 3 partition function. The duality predicts that the S 3 partition function is invariant under 
duality (up to an overall phase which is generally attributed to different framing [T2"|lo*5]). 
In our case, this statement translates to the following set of mathematical identities (for 

N < nk) 

i . / N \ N N 

U(Nt +1) ] = ^ / II II (2sinh(vrt l ,)) 2 II e^ + "*) = 



e 



53 L 

i#(N,k,n) Zss _ iV)^ 1 )] = (5.4) 

„MN,k,n) n fnk-N \ nk-N nk-N 



(nk — N)\ 




dtj] II (2sinh(vr^)) 2 n e e ^ +it 

i<j i,j=l 



where i?(iV, k, n) is some phase. In other words, the free energies F are duality invariant. 
We will test the validity of these identities first exactly in a simple, but non-topological, 
low-iV case, and then in the large-iV 't Hooft limit by using the results of our numerical 
computation. The matching in the latter case also provides supplemental evidence for the 
validity of these results and the accuracy of the numerical computation. 

Another way to derive the duality invariance of the free energy is by reduction of a 
superconformal index in four dimensions [21]. Starting from the electric-magnetic duality 
of an Sp(2N) theory in [M] it is possible to reduce the matching of the corresponding 
superconformal indices to the equality of three-dimensional partition functions in (15.41) , 

Other detailed aspects of the duality will be discussed in subsection 15.31 



5.1. A simple non-topological duality 



One can readily check the validity of the identities (15 .4p in the topological n — 1 case. The 
'hard' e^ ' -* factors drop out in this case and the integrals can be computed straightforwardly 
by brute force (see e.g. eq. (12.71) ). In general, one finds 
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U(k - Nt\ 



(5.5) 



We thank Grigory Vartanov for a discussion on this aspect. 
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Now let us consider a simple non-topological case. We set N = 2, k = l,n = 3 which 
leads to the duality 

U(2)¥ ] <+ U(l)^l . (5.6) 

The U(l) partition function involves a single Gaussian integral that gives (by analytic con- 
tinuation) 
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The U (2) partition function involves a double integral 
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We used the fact that e 
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With a simple change of integration variables we can 



factor out a Gaussian integral and reduce the computation to 
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(5.9) 



Using the identity sinh 2 (7rf) = cosh 2 (7rf) — 1 we can further split this integral into a Gaussian 
integral and an integral that can be computed directly with Mathematica 

1 7Ti 



Z s3 [U(2)l 
We conclude that 
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(5.10) 



(5.11) 



which verifies the prediction made by the duality in this case. 

Notice that the relative phase in (15. lip is not the same as that in eq. (15. 5p for k = 1. 
This implies that the relative phase has, in general, both a k and an n dependence. 

The particular example we have just discussed demonstrates the important role of the 
U(l) part of the gauge group in the duality. For an SU{2) theory at level 1 and no super- 
potential interaction a different type of duality was conjectured in 



5.2. Free energy matching tests at large-N 

We can provide additional tests of the above duality in the large- N 't Hooft limit by 
using the numerical results of section [2J Focusing on the planar part of the free energy we 
define 

f(X) = ^F(X t N) . (5.12) 
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Figure 5: Plots of the function D(\), defined in terms of the free energies in eq. (|5.14|) . for the A4 
theory (plot (a)), the A5 theory (plot (b)) and the Ag theory (plot (c)). The blue points represent 
the numerically computed results and the solid curves the duality-predicted functions X 2 (n — A)~ 2 . 



The duality predicts the equality 
F(X,N) = F I 
In Fig. we define 



D(X) 



& f(X)= Q-l) 2 /(n-A) . (5.13) 
f(n - A) 



/(A) 



(5.14) 



which according to the duality should be equal to \ 2 {n — A)" 2 . Using our numerical data 
we present how the matching works in three cases: n = 3, 4, 5 (i.e. the A 4 , A 5 , A 6 theories). 
For n = 3 the conformal window extends over the whole supersymmetric interval A G [0, 3]. 
For n = 4 the conformal window extends over the range A G (0.35,3.65). For n = 5 the 
conformal window extends over the range A G (0.65,4.35). In all cases we cover most of the 
conformal window fixing the probed range according to the ability of the numerical code 
to produce trustworthy results. The plotted data terminate at the self-dual point of each 
theory which lies at A s ^ = |. The agreement between the curve traced by the numerical 
points and the curve A 2 (n — A)~ 2 predicted by duality is rather convincing. 
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5.3. Other aspects of duality 

The successful matching of the free energies F in the above examples, combined with 
the arguments summarized in the beginning of this section, provide non-trivial evidence for 
the validity of the dualities (15.11) . It would be very interesting to show analytically the 
matching of the free energies f)5.4p for general N, k (modulo potential intricacies that are 
discussed below). Transformation properties of hyperbolic hypergeometric integrals [5oT455] 
may be very useful in this respect (as they have been in other cases, e.g. in the case of 
the Giveon-Kutasov duality [EHE2])- To the best of our knowledge the identities (j5.4p have 
not been proven in this way before (see [5H] for the appearance of analogous identities in 
four- dimensional theories). In anticipation of such results we would like to offer at this point 
a few additional comments on some particularly interesting, less discussed, physical aspects 
of the dualities (15. II) . For simplicity, we will focus on the large- N 't Hooft limit where the 
A n+1 theories are parameterized by n and the continuous 't Hooft coupling A. 

We have argued (see [22] for a detailed discussion) that there is a critical coupling A* +1 
where the deforming operator Tr X n+l becomes marginal. For A < A* +1 the deforming 
operator is irrelevant (for n > 3) and the theory erases the deformation in the deep IR where 
one recovers the undeformed A theory. For A > A* +1 the deforming operator is relevant and 
drives the theory to a new IR fixed point, which is the theory we are interested in. The 
global symmetry group of this theory is a single U(l), the R-symmetry group. The R-charge 
of the adjoint chiral superfield is fixed by the superpotential 

R(X) = -ij . (5.15) 

The duality is expected to match symmetries and operators. Indeed, this works very well 
inside the range of A's that we call the conformal window of this theory. This is the range 
where the deforming operator TrX n+1 is relevant both in the 'electric' U(N) theory and the 
'magnetic' U(nk — N) theory, namely when 

AG[A; +1 ,n-A; +1 ] . (5.16) 

In this case, the U(l) R-symmetry is identically the same on both sides of the duality because 
it is controlled by the same superpotential interaction, i.e. 

R(X)=R{X) = ^- (5.17) 
n + 1 

where X denotes the adjoint chiral superfield on the magnetic side. Notice, that because of 
this, the number of decoupled free chiral operators is the same on both sides of the duality. 
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Specifically, the single-trace operators 

TrX, TrX 2 , ■ ■ ■ ,Tr Xl^r\ (5.11 



n+l 



are free and decoupled on both sides of the duality ( [x] denotes the integer part of the number 
x). Accordingly, under duality not only F but also its modified version (14. 2|) is invariant. 

As an aside comment, we note the following interesting property of the RG flow from the 
UV undeformed A theory to the IR deformed theory A n+ i in this range of A. In the UV the 
number of free decoupled fields is m and this number is controlled by the exact R-symmetry 
of the A theory which is determined by F-maximization. Specifically, 



m 



2R(X) 



(5.19) 



Along the RG flow from A to A n+1 the R-charge increases and m drops in general from 
(I5.19P to ^-g- , which implies that a number of UV-free operators recouple in the IR and 
regain their interacting status. It would be useful to obtain a deeper understanding of the 
mechanism that realizes this effect. 

Returning to the A n+1 theory one can do a bit better in matching the operators. Following 
[6Ti] in appendix Owe demonstrate how one matches the general superpotential deformation 
( 15. 2p and the associated vacuum structure on both sides of the duality. A simple corrollary 
of that analysis is the map 

TrX* <* -TrX* , i = 1,2,. ..,n + l (5.20) 

for electric superpotential W = Tr X n+l and magnetic superpotential W = TrX n+1 . Under 
this map, a free decoupled operator from the list (I5.18P on the electric side maps to a free 
decoupled operator from the corresponding list on the magnetic side. 

This appears to be a satisfactory picture of how 3D Seiberg duality works in this class of 
Chern-Simons-Matter theories. Now one can ask what happens when A is increased further 
and dialed to lie in the window 

[n-\* n+1 ,n] (5.21) 

where a super symmetric vacuum continues to exist, but the deforming operator is relevant in 
the electric theory and irrelevant in the magnetic theory. F-maximization has allowed us to 
check that the size of this window is non-zero (this is essentially the statement that the exact 
R-charge curve, determined by F-maximization, lies below the upper bounding red curve in 
Fig.H]). This regime is analogous to the regime of the free magnetic phase in four-dimensional 
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M = 1 SQCD. The standard four-dimensional Seiberg duality is expected to work also in this 
phase by exchanging the 'very-strongly' coupled electric description to the IR free magnetic 
description. In our three-dimensional example we find a similar effect. In the vicinity of the 
supersymmetry breaking point, 0<n-A< 1, the superpotential deformation is so strongly 
coupled that it is in the verge of lifting the supersymmetric vacuum. In the dual description, 
where the dual 't Hooft coupling is small, A = n-A<l, the superpotential deformation is 
irrelevant. In this case the mismatching of symmetries, number of free decoupled operators 
etc. suggests that the correct R-symmetry cannot be read off the deforming superpotential 
in the electric description, and that one should rather use the magnetic description to study 
the theory in this regime (see [SO] for analogous comments in the case of the four- dimensional 
one-adjoint SQCD theory). 

It is interesting to understand how this story implements itself on the level of the free 
energies F computed via localization in the regime (15.211) . A naive possibility that seems 
unlikely to work (for the above reasons) is to try to match the free energy computed via 
localization at R = in the electric theory and the free energy computed via localization at 
R determined by F-maximization in the A theory on the magnetic side. Another possibility 
is that the mathematical identities that give rise to (15. 4p are blind to the existence of a 
conformal window and apply for all N, k in the supersymmetric interval. In that case, 
outside the conformal window the matrix model F does not refer to the actual free energy 
of the theory. Yet another possibility is that they only apply for N, k within the conformal 
window, but then it is interesting mathematically to identify precisely how they break down 
outside of it. 

6. Tests of the F-theorem at large-A 

Recently, an analog of the c-theorem in two dimensions has been put forward for three- 
dimensional quantum field theories in Refs. [H1I21]- The precise proposal, which is known 
as the F-theorem, states that F (the free energy of the S 3 partition function we have been 
computing) always decreases along RG flows. In the absence of a general proof it is instructive 
to extend the tests in as many classes of theories as possible. In this section we test the 
validity of the theorem in a new set of RG flows that occur within the setup of the A and 
A n+1 theories. 
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6.1. A web of RG flows 

By adding relevant superpotential interactions to the one-adjoint A CSM theory one can 
generate a web of RG flows. We have already discussed some of these flows. In the range 
of coupling where a chiral ring operator TrX n+1 is relevant (there is such a range for any 
n > 1), one can use the operator to deform the A theory by the corresponding superpotential 
and generate a supersymmetric RG flow towards a new IR fixed point that we called A n+1 . 

Another general class of RG flows arises by using an A n+ i theory as the UV fixed point 
and deforming it by a general polynomial superpotential of the form 

n a- 

i>% TV- yn+l—i 



W = Y - TrX n+1 -* . (6.1) 



i=0 

The supersymmetric vacua of this theory can be found by solving the F-term equations 



n 



w'(x) = Y,g l x n - l = goIl( x -^P > E^ = w ( 6 - 2 ) 

i=0 i=l i=l 

where dj are parameters directly related to the coefficients gi of the superpotential poly- 
nomial. In the far infrared the theory flows to a product of A nj+ i theories specified by a 
partition of the total rank iV 

f/(AT)i" +1) U^)^ x • • ■ x U{N t )t l+l) , J2 N i = N ■ ( 6 - 3 ) 

i=l 

There is a supersymmetric vacuum for each of the factors in this product provided Ni < riik 
for alH = 1,2,..., £. 

6.2. General predictions of the conjectured F -theorem 

For each of the above RG flows the F-theorem predicts the inequality 

F uv > F IR . (6.4) 

For the general flow f)6.3p this inequality implies 

^ fe +1) >E4:: fe +1) - (6.5) 
i=i 

Focusing on the planar contributions /(A; n) in the large- iV limit (see eq. (15.121) ) we obtain 

/fe >E^M;^) ( 6 - 6 ) 
\i=i i=i / i=i 
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where by definition 

t 

N l=Xl N , \ l = x l \, 0<Xi<l, j>i = 1 . (6.7) 

i=l 

Notice that in these flows it is in general possible to start in the UV from an A n+1 theory 
inside its conformal window and end with a product of theories where some of them lie 
outside their conformal window. In these theories the deforming operator may be irrelevant 
in which case the true IR fixed point is the A theory and the free energy should be computed 
accordingly. 



6.3. Tests of the F-theorem 

We proceed to demonstrate the validity of the above inequalities in a set of examples 
with increasing complexity. We will concentrate on the large-iV limit where we can make 
use of the numerical results of section [2J 



6.3.1. A ->■ A n+l! A n+1 -»■ A^z+i, ri < n 

We begin with the RG flow A — > A n +i. We pick A < n to be such that the deforming 
operator Tr X n+1 is relevant in the UV theory. In these flows 

2 

Ruv < Rir = ■ (6.8) 

n + 1 

Since F(Ruy) (determined by F-maximization) is a global maximum in all these cases, the 
inequality (16.41) immediately follows. 

As a simple example of the general flow (16. 3p we can consider the case 

U{N)l uv+1 U(N) n k IR+1 , n IR <n uv , N< n IR k . (6.9) 

In this case, the 't Hooft coupling A is constant along the RG flow. We choose it so that it 
lies inside the conformal window of the UV theory, namely A G [Xn uv+ i, nuv — Ki uv +i]- Since 
^n IR +i < Kiu +i: ^ foll° ws that t ne IR theory will also be within its conformal window. As 
a result, 

Ruv = 2 — < Rir = 2 — ■ (6.10) 

nuv + 1 n IR + l 

The monotonicity of F as a function of R (at fixed A) guarantees again the F-theorem 
inequality (16.41) . 
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6.3.2. A 8 ->A 4 (g>A5 

We move on to a more intricate one-parameter family of examples 

(n = 7, A) -»■ (n = 3, Ax) ® (n = 4, A(l - x)) , 0<x<l. (6.11) 

Further constraints on the parameter x may arise by requiring that the IR theories have a 
supersymmetric vacuum. Being obvious we will not write out these constraints explicitly 
below (in the specific cases that we consider they are in any case vacuous). The range of the 
conformal windows of the A 4 , A 5 , A 8 theories are 

A 4 : (0,3) , A 5 : (0.35,3.65) , A 8 : (1.15,5.85). (6.12) 

For illustration purposes let us take A = 2 (similar results can be obtained also for other 
values of A). The F-theorem inequality (16. 6p becomes in this case 

A(x) = /(2; 7) - x 2 f(2x; 3) - (1 - x) 2 f(2(l - x); 4) > . (6.13) 

The results of the numerical computation are depicted in plot (a) of Fig. EJ We observe that 
the inequality (I6.13P is verified. 

6.3.3. A 10 ^A 6 <g>A 5 

The results of a similar computation with different values of n are depicted in plot (b) of 
Fig. [6J Now we consider the one-parameter family of flows 

(n = 9, A) -» (n = 5, Ax) <g> (n = 4,A(1 - x)) , 0<x<l (6.14) 

and for concreteness we set A = 3. The range of the conformal windows of the A 5 , A 6 , A 10 
theories are 

A 5 : (0,35,3.65) , A 6 : (0.65,4.35) , A 10 : (1.7,7.3). (6.15) 

The F-theorem inequality (16.61) becomes in this case 

A(x) = /(3; 9) - x 2 /(3x; 5) - (1 - x) 2 /(3(l - x); 4) > . (6.16) 

Again, it is apparent from the results depicted in plot (b) of Fig. [6] that the inequality (16.191) 
is verified. 
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Figure 6: P/ois expressing the function A m i/ie case of the three RG flows described respectively 
in subsections 1 6. 3. ffi 1 6'. ff. iJ[ and \6. Sl\ The F-theorem requires that A > 0, which is indeed verified. 
The results are based on the numerical solutions of section^ 

6.34. A n -> A 4 <S> A 4 ® A 5 

Finally, we report a test of the F-theorem for a two-parameter family of RG flows defined 

as 

(n = 10, A) ->■ (n = 3, A21) ® (n = 3, Ax 2 ) <8> (n = 4, A(l - 21 - x 2 )) (6.17) 

with < Xi t 2 < 1, < 1 — Xi — x 2 < 1. In the case depicted in Fig. ^c) A = 3. The range 
of the conformal windows of the A 4 , An theories are 

A 4 : (0,3) , A n : (2,8) . (6.18) 

The F-theorem inequality (16.61) is more complicated in this case 

A(xi,x 2 ) = f(3;10)-x 2 1 f(3x 1 ;3)-x 2 2 f(3x 2 ;3)-(l-x 1 -x 2 ) 2 f(3(l-x 1 -x 2 );A) >0. (6.19) 

The results depicted in plot (c) of Fig. verify this inequality. 
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One could proceed in this way to generate more tests of the F-theorem. 
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Appendices 



A. Multi-cut saddle point configurations 



Our purpose in this appendix is to demonstrate that the saddle point equations (12.41) 
have general multi-cut solutions. We will discuss the presence of these solutions in a pertur- 
bative expansion of the equations in A 1 focusing, for illustration purposes, only on the 
leading term of the expansion. Our main concern here is not to find explicit solutions to 
the perturbative equations, but to demonstrate the allowed types of expansions. We have 
verified the existence of explicit solutions numerically. 

A perturbative analysis of the saddle point equations 

(1 - R) sinh(27rt ii ) + t y sin(27r.R)' 



Z 



-u 



A 



N 



E 



coth(7rti 



cosh(27r£jj) — cos(2ttR) 
together with the R-extremization equation 



0, 



1,2, 



Re 



N 

E(i 



R + itij) cot 7r(l — R + itij) 







,N 
(A.l) 

(A.2) 



appeared in [5] . In this appendix we are interested in general values of R, hence we will not 
worry about the R-extremization equation. 
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In order to keep the discussion simple and compact we begin by listing a few characteristic 
expansions. 



The standard one- cut solution crossing the origin 

This expansion is based on the ansatz 

U = VX(4 0) + a4 1) + ---) , i = 1,2,... ,7V. (A.3) 

Inserting this expansion into the saddle point equations (1A.1I) and expanding in powers of A 
we get 

+ . . . = . (A.4) 

The dots represent higher powers in the expansion. Ref. [5] explains how one can proceed to 
solve these equations systematically at any order. At leading order the equations (which are 
independent of R) reduce to those of the familiar Wigner model and the eigenvalues exhibit 
a semicircle distribution. 

This expansion captures the one-cut solutions that reproduce correctly the perturbative 
field theory computation of the R-charge and are the key players in the numerical analysis 
of the main text. 



1 



if 



(o) 



N 

E 



7tN t (0) 



Two-cut solutions crossing the imaginary axis at ± 
Consider an expansion based on the ansatz 



Of 1 



U = ~ + A4 0) + A 2 *?' , i = ^ + 1,^ + 2 NeS~. (A.5) 

Inserting this expansion into the saddle point equations we find for m G IR + — N 



and for m G N 



1 / m 1 A 



I - ± = A^2 + ^gfF^wJ+'"= ' < A7 » 

One can show that only the m G N case admits solutions. We have not attempted to 
find these equations analytically, but have verified numerically that they exist (see e.g. plot 

(b) of Fig. up. 
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Two-cut solutions crossing the imaginary axis at -j- 
In this case we consider the ansatz 

(1= CH + Hi + At m + AH m i i=1 , 2 ,...,^ €S+ , 

ij = _(fl^)l +At <., +AWi j = £ + 1 ,£ + 2 ,..., WeS -. (A.8) 

Inserting this expansion into the saddle point equations we obtain to leading order 

_ 1 / R + m 1 x , 1 m + 1 ^ 1 \ . 

General multi-cut expansions 

One can also entertain the possibility of more general multi-cut solutions based on the 
following ansatz 

N eigenvalues in So : U = \ a tf + . . . , 

iV m eigenvalues in <S* : t« = ±— — h At- -* + . . . , (A. 10) 

M m eigenvalues in S±' R : ^ = ± ( jR + m ) 1 + ^(o) + . . . ? 



so that 



N + J2 2 ( N m + M m ) = N . (A.ll) 



m=l 



a = ~, 1 are exponents with special features but other values may also be allowed. 

We have not attempted to solve the above equations analytically, but have verified nu- 
merically the existence of multi-cut solutions of the above more general type. In general, 
the actual solution of these equations will impose further constraints on the allowed values 
of the integers No, N m , M m . 

The above list covers the full set of allowed expansions. 



B. Spontaneous SUSY breaking implies Z$ 3 ' = 0: outline of an argument 

In this appendix we outline a potential line of reasoning in favor of statement (6) in 
section Our purpose is to show that in M = 2 theories with spontaneous supersymmetry 
breaking the Q-deformed partition function Zg° c ^ vanishes. We describe on general grounds 
why this is a potentially correct generic property and highlight some of the key aspects that 
need to be clarified in order to reach a proof. 
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Let us begin with a few general comments. For a theory on a 3-manifold M. with an S 2 
boundary there is a natural Hilbert space % associated with the boundary S 2 [SB] • The path 
integral of the theory over M. with boundary conditions $|g2 = \ 

Z(x)=f e~ s (B.l) 

"'*| S 2=X 

defines a vector \x) in the Hilbert space %. We are interested in the case where M. is a 
semisphere. 

With the insertion of a generic operator O at the pole of the semisphere the path integral 
with boundary conditions \ gives the amplitude (0\x)- This formalism is the basis of the 
operator-state correspondence (see, for example, [SI] for a review). In a bra-ket language 

Z( X ) = (m (B.2) 

where \Q) is the vacuum state of the theory corresponding to the insertion of the identity 
operator. 

The partition function on S 3 can be written in the following form 



z sa = (fi|n> = / d x z{x)z\ x ) = J d x (n\ x ){x\n) 



(B.3) 



The second equality corresponds to gluing two 3-spheres at the S 2 equator and summing 
over the boundary conditions x- We will be applying (1B.3|) to the Q-deformed path integral 
(13. ip . In this language 

dZ s s(t) 



■<n| J{Q,v}\n) . (b.4) 



dt 

As a result, in a theory with spontaneously broken supersymmetry Q\Q) ^ and generically 
the derivative dZs ^ t — is non-zero. From now on we will be referring exclusively to the t = oo 
theory whose partition function is denoted as Z^ c ^ in the main text. 

Another useful piece of information that we need is the following. Localization is based 
on a supercharge Q. There is an additional supercharge with anticommutator 

{Q,Q^} = M + R (B.5) 

which is also a symmetry of the theory. M, which is part of the USp(2, 2) conformal group 
on S 3 , is a rotation on S 3 that can be viewed as a translation along the Hopf fiber of 
S* 1 S 3 — > S 2 . R is the R-symmetry operator, 



15 



For additional details we refer the reader to Ref . [5] . 
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The ground states f2 have zero R-charge. Hence, using a standard argument based on 
the correlation function 



(Q\M\Q) = (Q\{Q,Q^}\Q) = \Q\Q)\ 2 + Q^\Q) 2 (B.6) 

we deduce that a ground state is supersymmetric if and only if it has M = 0. Notice, however, 
that the second equation in flB.61) is strictly valid in Lorentzian signature. In Euclidean 
signature is not the Hermitian conjugate of Q, but an independent supercharge. In what 
follows we will nevertheless assume that the above conclusion, i.e. that a ground state is 
supersymmetric if and only if it has M = 0, is also correct in the Euclidean theory. Clearly, 
this is a point that deserves further justification. 

Furthermore, one can show, using an argument in [35], that the zero-M ground states 
are precisely the zero R-charge states \a) that obey the equation 

Q\a) = 0, \a)^Q\(3) (B.7) 

for any state This concludes our short parenthesis of introductory comments. 
By definition, Zg° c ^ receives contributions from states \x) that solve the equation^ 

{Q,V}\ X ) = 0. (B.8) 

Applying Q to this equation and using Q 2 = we find that the states VQ\x) are annihilated 
by Q, i.e. that 

QVQ\x) = . (B.9) 

Therefore, according to the above general discussion, if the states VQ\x) are zero R-charge 
states that are not Q-exact we may infer that they have M = and are therefore supersym- 
metric ground states. Then (b) follows naturally. When supersymmetry is spontaneously 
broken there are no supersymmetric ground states. Hence, there are no zero R-charge, zero- 
M ground states, including states of the form VQ\x), and therefore no states x that obey 
equation (1B.8|) that can contribute to the path integral Zg° c \ As a result, Zg° c ^ receives 



no contributions from the sum / dx in the rightmost expression in eq. (IB. 31) and vanishes 
identically as stated in {b). 

The assumption that the states VQ\x) have zero R-charge is plausible for the following 
reason. The implementation of localization in general Yang-Mills theories with arbitrary 



16 More specifically, since localization forces the path integral on configurations $ with the property 
{Q, y}$ = the same property must be obeyed by continuity by the boundary conditions x- 
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matter [IH3] and the eventual reduction of the path integral to a matrix integral over the 
real, zero R-charge scalar field a in the M = 2 vector multiplet, suggests that, at least in 
this general class of theories, x is a state with zero R-charge. Since the R-charge of VQ is 
zerd_j we conclude that VQ\x) is a zero R-charge state as well. In more general situations 
one has to examine precisely how localization works. 

The other part of the argument requires showing that the state VQ\x) is not Q-exact, 
namely that there is no state \j3) for which 

VQ\x) = Q\P) ■ (B.10) 

This property seems plausible, but we have not been able to find a rigorous proof. 

It would be interesting to know if the specific choice of V plays any particular role in this 
argument. In general, Zg° c ^ is expected to depend on the choice of V when a supersymmetric 
vacuum is absent, but the property Zg° c ^ = may not. The latter would have to be true 
for (b) to hold in its current form, otherwise one would have to specify a special class of V 
functionals for which (b) is valid. 

Moreover, part of our proposal is that any unknown trial R-charges in Zg° should 
be fixed by l-Z^ ^ -minimization even in the absence of a supersymmetric vacuum. This 
prescription is motivated by the requirement to have a universally-prescribed quantity that 
reproduces the physical sphere partition function when supersymmetry is not broken. The 
implications and necessity of this prescription for the validity of (b) should be clarified. 



C. Matching deformations in the A n+ i duality 

The duality (15. ip is expected to hold for arbitrary superpotential deformations of the 
type (16.1 1) . Requiring the general matching of the vacuum structure on both sides of the 
duality gives useful information about the precise map between chiral ring operators. We 
proceed to derive this map by suitably adapting the four- dimensional analysis of [BP] . 

We define the 'electric' theory as the U(N)^ +1 ^ theory deformed by the general super- 
potential 

n n 

W = Y Tr X n+1 ~ l . (C. 1) 

The dual 'magnetic' theory is a U(nk — N)^ 1 ^ theory deformed by a dual superpotential 



17 Indeed, the R-charge of the Q-deformation in (|3.ip must be zero. The explicit formulae for V in [IH3] 
verify this property. 
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of the form 

n a- 

W = T, J: ■ Tr Xn+l ~ l + <*({&}) ■ (C2) 

X denotes the chiral superfield in the adjoint representation of the dual theory. We have 
allowed for a constant term in the superpotential denoted as a({g,i\). Our purpose is to 
determine g~i and a as functions of the electric superpotential couplings gi. 



Electric theory 

Consider, for example, the case where 

n n 

W\x) = E^""' = 9o II - a i) (C3) 

i=0 i=l 

with all dj different. As described in the main text, in the deep IR the theory splits into a 
set of decoupled A 2 theories with gauge groups U(Nj) labelled by the sequence of integers 

n 

N 1 <N 2 <---<N n , Y, N * = N - (C- 4 ) 

i=l 

Ni eigenvalues of the matrix X reside in the i-th minimum of the potential V = \W'(x)\ 2 
labeled by Oj. Each of these A 2 theories has a supersymmetric vacuum if and only if Ni < k. 



Magnetic theory 

Similarly, on the magnetic side 

n n 

W'{x) = Y^gix"-* = g U(x - (C.5) 

i=0 i=l 

with all a, different. In the deep IR the theory splits into n copies of the A 2 theory with 
gauge group U(Ni), such that 

n 

Ni < N 2 < ■ ■ ■ < N n , ^TNi = N = nk — N. (C.6) 

i=i 



The duality must work individually for each of the n components of the IR theory, hence 

Ni = k-Ni, z = l,2,...,n. (C.7) 

For this property to be realized with any choice of partitions {N} the electric and magnetic 
superpotentials must be closely related. In particular, whenever any number of critical 
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points a, coincides, the same number of dual critical points must coincide as well. This 
requirement puts constraints on the dual couplings g^. 
A potential solution to these constraints is to set 

9i = cgi , % = 1,2, . . .,n (C.8) 

where c is some constant. We adopt this solution and, following [60J, fix our conventions so 
that c = — 1. 

The duality map gx — > gi is closely related to the map between the chiral ring operators 
TrX* and TrX J . More information about this map can be inferred in the following way. 
Define the free energy F(^) of the theory in flat space as 

e - / d A xd 2 e¥{ ai )+c.c. _ / - / d 3 xd 2 eW(X, gi )+c.c.\ g\ 

where gi are background chiral superfields. Then 

1 r)F 

(TrX n+1 ^) = ^. (CIO) 



n + 1 — % dg. 
Similarly, on the magnetic side 

1 - BF 

— : (Tr X n+1 ~*) = — . (Gil) 

n + 1 — i ogi 

The duality implies 

F(g t (g)) = F( 9l ) . (C.12) 

Hence, 

TrX n+1 ~ i = y TrX n+1 -i$± + ^ . (G13) 

n + l — i j^n + l-j dgi dgi 

Taking the vev of both sides of this equation we find additional constraints on the map 
9i —I 9i- Most notably, the vevs of the left hand and right hand sides, which depend 
non-trivially on the particular vacuum (namely the partition {iVj}), must satisfy a relation 
independent of the particular vacuum. Implementing the ansatz (1C8|) (with c = —1) we 
obtain the relation 

— 

TlX n+i-i = _ TiX n+i-i + ^ n + 1 _ •) i = l,2,...,n. (C.14) 

@9i 

Taking the vev of both sides of this equation in a vacuum specified by the partition {iVj} 
and the critical points aj = a, gives (taking also the duality eqs. (1G7I) into account) 

ti -\- 1 — z d(y 

w n+ i_i = — (G15) 

k dg { 
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where following [BU! we have defined 



n 

", -»<)•' • (C16) 



Solving ( 1C.15I) one finds 



n + 1 ~[ V n + 1 — ij 



The following identities are useful in deriving this result 

d ( u i+m \ _ . . ™ { duj 



independent of i , V igi^- = juj . (C.18) 
dg { \i + mj fri \ dgij 

We conclude that the operator map takes the form 

TiX n+1 ~ l = -TrX^^ + kUn+^i , i = 0,l,...,n (C.19) 

and the superpotentials W, W match with the identifications (1C.8I) and ( 1C.17I) . In the 
special case where gi = (i — 1, . . . , n), go ^ the map flC.19j) becomes simply 



= _TrX" +1 -\ i = 0,l,...,n. (C.20) 
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